Photon interferometry of Au+Au collisions at the Relativistic Heavy-Ion Collider 
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We calculate the two-body correlation function of direct photons produced in central Au+Au 
collisions at the Relativistic Heavy Ion Collider. Our calculation includes contributions from the 
early pre-equilibrium phase in which photons arc produced via hard parton scatterings as well as 
radiation of photons from a thermalized Quarli-Gluon-Plasma (QGP) and the subsequent expanding 
hadron gas. We find that high energy photon interferometry provides a faithful probe of the details 
of the space-time evolution and of the early reaction stages of the system. 



Recent data from the Relativistic Heavy-Ion Collider 
(RHIC) at Brookhaven National Laboratory have pro- 
vided strong evidence for the existence of a transient 
QGP - among the most exciting findings are a strong 
collective flow saturating the hydrodynamic limit 0, , 
a suppression of high-p^ hadrons indicative of parton en- 
ergy loss i^J , and an evidence for hadronization by re- 
combination of quarks from a deconfined phase . Being 
of hadronic nature, these probes are necessarily affected 
by the late stages of the reaction and do not provide 
unambiguous information about the early equilibration 
phase. Photons and dileptons, on the other hand, de- 
couple from the system immediately after their creation, 
and thus serve as complementary deep probes of the hot 
and dense matter. Recent studies have suggested that 
rescattering of partons enhances the production of pho- 
tons HQ. 

One problem affecting photons is that the inclusive 
photon spectrum is dominated by radiative decays of 
long-lived hadrons (especially tt" and rj). Although the 
subtraction of these sources is possible with great effort, 
this problem severely limits the use of photons as pene- 
trating probes. Two-photon momentum correlations do 
not suffer from this problem, because the geometric size 
of the distribution of decay vertices is extremely large 
on a nuclear length scale, and thus all correlations be- 
tween photons from the decay of two different hadrons 
are confined to extremely small relative momenta. This 
property makes photon interferometry a promising tool 
for the experimental study of the pre-equilibrium dynam- 
ics of relativistic heavy ion collisions. 

In the present work we employ the parton cascade 
model (PCM) [3,01 and ideal hydrodynamics to cal- 
culate two-body quantum correlations of high en- 
ergy photons and show that they can provide a sensi- 
tive test of the space-time evolution of the matter in 
the collision zone. The peculiar features of the corre- 
lation function make it possible to differentiate between 
direct photons emitted in the early pre-equilibrium re- 
action phase from those produced at later times in a 
thermalized QGP and subsequent hadron gas. Follow- 
ing the early suggestion of Neuhauser photon in- 



terferometry of the QGP was investigated extensively by 
several authors 0, Q, 0, 0, 0, 0| . However those 
studies focused only on photons emitted from a thermal- 
ized source. Recently, interferometry measurements with 
photons of transverse momenta of 100-300 MeV/c have 
been reported at SPS energies [l^ . 

Here we are primarily interested in the higher- 
energy photons originating in the rapidly evolving pre- 
equilibrium phase when the initial-state partons undergo 
multiple inelastic collisions. We use the Parton Cascade 
Model to describe this phase of the collision. The basic 
assumption underlying the PCM is that the state of the 
dense partonic system can be characterized by a set of 
one-body distribution functions Fi{x^ ,p°'), where i de- 
notes the flavor index {i — g, u, u, d,d, . . .) and a;^,p" are 
coordinates in the eight-dimensional phase space. The 
time-evolution of the parton distribution is described by 
a relativistic Boltzmann equation, whose collision term 
includes all scattering processes involving quarks, glu- 
ons, and photons in lowest-order quantum chromo- and 
electro-dynamics. A low momentum-transfer cut-off p™'" 
is introduced to regularize the infrared divergence of the 
perturbative parton-parton cross sections. Additionally, 
we include the branchings q — > qg, q — > 57, g gg 
and g qq following two-body collisions in the leading 
logarithmic approximation. The soft and collinear singu- 
larities in the showers are eliminated by terminating the 
branchings when the virtuality of the time-like partons 
drops below /io = 1 GeV. The present work is based on 
the thoroughly revised, corrected, and tested implemen- 
tation VNI/BMS of the parton cascade model We 
already presented the rapidity distribution and the trans- 
verse momentum spectra of direct photons predicted by 
VNI/BMS in an earlier publication 0. Here we explore 
two-photon correlations for the same model. 

The PCM calculation of the hard pre-equilibrium pho- 
tons is augmented by a hydrodynamical calculation for 
photon emission from a thermalized QGP and subse- 
quently forming hadron gas. Here, we employ a boost 
invariant hydrodynamics [iflj, which has been used ex- 
tensively to explore photon and dilepton production and 
also hadronic spectra. The initial conditions are esti- 
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FIG. 1: Upper panel: The production rate (per event) of 
hard photons in a central collision of gold nuclei, at -^sjvjv ~ 
200 GeV as a function of time in the center-of-mass system. 
Lower panel: Spectrum of photons from various sources: QM 
denotes photons emitted from a QGP and HM denotes pho- 
tons emitted from a hadron gas. 



mated from a particle rapidity density of about 1260 near 
y = for a central collision of gold nuclei at y^sj^ = 200 
GeV. We assume that a chemically and thermally equi- 
librated plasma is formed at tq = 0.3 fm/c. We use the 
complete leading-order results for the production of pho- 
tons from the QGP 0| and the latest results for the 
radiation of photons from a hot hadronic gas |2l| . 

The upper frame of Figure ^ shows the emission rate 
dN^/dt for photons with a transverse momentum of 
(1 ± 0.1) GeV/c around mid-rapidity as a function of 
time, i = is defined in the nucleus-nucleus center-of- 
mass frame, as the instant when the two Lorentz con- 
tracted nuclei completely overlap. 88% of all photons 
produced via hard parton scattering in the PCM are 
emitted within the first 0.3 fm/c of the reaction. Be- 
yond this time, due to its artificial infrared cutoffs, the 
PCM becomes an increasingly ineffective description of 
the medium evolution and does not achieve or maintain 
full equilibrium among the interacting partons. Assum- 
ing an initial time on the order of 0.3 fm/c, the hydro- 
dynamic calculation allows for a smooth continuation of 



the emission rate as a function of time for the QGP and 
subsequent hadron gas phase. Note that the production 
of photons due to the primary-primary scatterings only 
is approximately symmetric around t — 0, confirming the 
build-up of secondary scatterings over time. 

The relative weight of photons emitted in the pre- 
equilibrium phase from the PCM vs. those stemming 
from a thermalized QGP and subsequent hadron gas can 
be studied as a function of transverse momentum in the 
lower frame of Figure ^ For a transverse momentum 
of fcr = 1 GeV/c pre-equilibrium and thermal photons 
both contribute roughly 50% to the total photon yield. 
We therefore expect both contributions to be of impor- 
tance for the calculation of the correlation function. At 
kx — 2 GeV/c, however, direct photons from the PCM 
outweigh thermal photons by one order of magnitude. 
Here we expect only negligible contributions of thermal 
photons to the correlation function and a clear signal of 
the early photon emitting source. 

The correlation between two photons with momenta 
ki and k2, averaged over their spins is 



where S{x, k) is the photon source function for a com- 
pletely chaotic source, and 

q = ki-k2, K = (ki+k2)/2 . (2) 

Our assumption of a chaotic source neglects contribu- 
tions to the 2 photon yield by correlated 2 photon emis- 
sion. We estimate these contributions to be small and 
not to influence the shape of the correlation function. 

In the parton cascade model calculation, where the 
probability distribution of the production vertices and 
momenta of the photons are known, the correlation func- 
tion C is calculated by following the Wigner function 
scheme of [l^ using the Hansa code ^23], i.e. by rewrit- 
ing the numerator of the second term in Eq.(^ as, 

Num(q,K) =^cos[q-(xi-Xj)]. (3) 

where the sum runs over the pairs whose individual mo- 
menta lie in a small bin around K of width e in each of 
the four directions. 

The emission vertices (x, k) generated by a semiclassi- 
cal transport model or the monte-carlo sampling of emis- 
sion rates generally do not represent a valid Wigner func- 
tion, since the uncertainty relation Ax Ak > ?i is not in- 
corporated. In order to remove this deficiency we have 
smeared the emission vertices by add ing random fractions 
of h/pi) to the emission vertices Xi |2J|. 

We shall present our results for the correlation function 
C(q, K) in terms of the so-called outward, sideward, and 
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FIG. 2: The outward, sideward, and longitudinal correlations 
of direct photons predicted by the parton cascade model, at 
Kt ~ 2 GeV/c. Note the change of scales along the x and y- 
axes for the longitudinal correlation. The inclusion of thermal 
photons changes the results only marginally. 



longitudinal momentum differences. We represent the 
photon four-momentum in the form 

fc'* — {kx cosh y, kx cos?/;, fcr sini/;, kx sinh y) (4) 

where kx is the transverse momentum, y is the rapidity, 
and ?A is the azimuthal angle. Using the same notation 
for the relative and the average momenta we can 
write (251] : 

giong = \kiz- k2z\ = IfciTsinhyi - fc2Tsinhy2| (5) 
(Jout = qx • Kx/i^T (6) 
gsidc = IQT - qoutKT/i^Tl . (7) 

The corresponding radii for a completely chaotic source 
are often obtained by fitting the correlation function 
C(gout, ^sidc, qiong) to the functional form: 

l + (l/2)-exp [- (q^^t^out + 9sidc-Rsidc + 'Zlong^ng) /2 ■ 

Note that this provides that (qf) = where i=out, 

side, long. 

Transport calculations have employed several proce- 
dures to improve the significance of their results in order 



to compare with experimental data. Either the source 
function is parametrized, permitting an analytical eval- 
uation of eq. (Q), or the particles generated in the cal- 
culation are "dispersed" over many particles to enhance 
the pair statistics. For our analysis we have accumulated 
125000 events from VNI/BMS to make an accurate deter- 
mination of the source radii possible by direct evaluation 
of cq. ©. We supplement these results with photons 
from a Monte-Carlo sampling of the hydrodynamics cal- 
culation to obtain the production vertices and momenta 
of the photons, which are combined with those from the 
PCM in our correlation analysis. 

In Fig. 121 we show our results for photon momenta of 
kx — (2±0.2) GeV, where the high energy photons stem- 
ming from hard parton scattering completely dominate 
the total photon yield. We study the three-dimensional 
correlation function, in terms of one momentum coordi- 
nate, limiting the other two momentum coordinates to 
zero. The sideward radius as inferred from photons is 
^sido ~ 3.8 fm, while the corresponding outward radius 
is found to be i?out ~ 3.7 fm. These values are some- 
what smaller than the Gaussian equivalent radius of the 
density distribution of a Au nucleus, refiecting the ge- 
ometric shape of the colliding nuclei. The very small 
difference between the two radius parameters indicates a 
short source lifetime. 

The longitudinal correlation radius is much smaller, 
^long ~ 1-6 fm. This is not unexpected 0| and has its 
origin in the fact that these photons are mostly emitted 
when the longitudinal extension of the system is quite 
small and the velocity gradients are large. We have ver- 
ified that C(giong) is robust with respect to reasonable 
variations of the center of the rapidity window and its 
width, which determines e. A measurement of this value 
of i?iong would clearly identify the source of these photons 
as prehadronic. 

In figure |3| we show our results for photons at kx = 
(1 ± 0.1) GeV. For this choice of kx pre-equilibrium pho- 
tons and thermal ones from the QGP and hadron gas 
phases of the reaction contribute about equally to the 
yield - therefore both contributions need to be taken into 
account when calculating the correlation functions. As 
is to be expected, the radii extracted from the correla- 
tion functions lie between those of the individual contri- 
butions (thermal photons vs. pre-equilibrium photons), 
however, they in general do not reflect the average of the 
radii as the relative single photon yields would suggest, 
due to the influence of the cross-term, i.e. the interfer- 
ence of thermal photons with pre-equilibrium photons. 

For low kx, thermal photons should completely domi- 
nate: at kx ~ 0.2 GeV, our calculations predict i?sido of 
about 5 fm, with i?out/.Rside ~ 2. 

In summary, we have calculated the correlation func- 
tions of high energy photons emanating from the pre- 
equilibrium stage as well as the thermal QGP and hadron 
gas stages of the matter formed in central Au+Au colli- 
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FIG. 3: Outward, sideward and longward intensity correla- 
tion of photons at 1 GeV, considering only PCM (BMS), only 
thermal (hyd), and PCM+thermal photons (BMS+hyd). 



sions at RHIC. Our study predicts that photon intcrfer- 
ometry will reveal a small source of brief duration for pho- 
tons at transverse momenta kx > 2 GeV/c due to pre- 
equilibrium emission. This contrasts with a much larger 
source, revealed by thermal photons at lower momenta, 
as they mostly originate from the QGP and hadron gas 
stage when the system has already expanded signifi- 
cantly. 
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